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Embedded Si and Ge nanocrystals (NCs) in wide band-gap matrices are studied theoretically using 
an atomistic pseudopotential approach. From small clusters to large NCs containing on the order 
of several thousand atoms are considered. EfTective band-gap values as a function of NC diameter 
reproduce very well the available experimental and theoretical data. It is observed that the highest 
occupied molecular orbital for both Si and Ge NCs and the lowest unoccupied molecular orbital for 
Si NCs display oscillations with respect to size among the different irreducible representations of 
the C-iv point group to which these spherical NCs belong. Based on this electronic structure, first 
the interband absorption is thoroughly studied which shows the importance of surface polarization 
effects that significantly reduce the absorption when included. This reduction is found to increase 
with decreasing NC size or with increasing permittivity mismatch between the NC core and the 
host matrix. Reasonable agreement is observed with the experimental absorption spectra where 
available. The deformation of spherical NCs into prolate or oblate ellipsoids are seen to introduce 
no pronounced effects for the absorption spectra. Next, intraconduction and intravalence band 
absorption coefficients are obtained in the wavelength range from far-infrared to visible region. 
These results can be valuable for the infrared photodetection prospects of these NC arrays. Finally, 
excited-state absorption at three different optical pump wavelengths, 532 nm, 355 nm and 266 nm 
are studied for 3- and 4 nm-diameter NCs. This reveals strong absorption windows in the case 
of holes and a broad spectrum in the case of electrons which can especially be relevant for the 
discussions on achieving gain in these structures. 

PACS numbers: 73.22.-f, 78.67.Bf, 78.40.-q 



I. INTRODUCTION 

The field of silicon and germanium nanocrystals (NCs) 
is very active due to important technological achieve- 
ments and prospects particularly in connection with op- 
tics such as light emitting diodes and lasers fii^i^ Two 
fundamental processes describing the interaction of light 
with matter are the photon absorption and emission. 
In the context of NCs, it has been shown both exper- 
imentally^ and theoretically^"^ that the interface prop- 
erties have dramatic effects on the emission properties. 
On the other hand, absorption measurements are less 
sensitive to surface quality and allow for a more di- 
rect characterization of the intrinsic structure of NCs^- 
Therefore, the study of the direct photon absorption in 
NCs can provide a clear physical understanding. More- 
over, with the ever-growing importance of renewable en- 
ergy resources, the research on the new-generation pho- 
tovoltaics has gained momentum and hence the sub- 
ject of direct photon absorption in nanocrystalline sili- 
con (nc-Si)i^ However, both experimentall^ r^^i^^'^^i^^'^^^i^^ 
and theoretically )^^i^''i^^i^^i^°i^^i^^i^^ researchers till now 
have predominantly focused on the interband absorption 
process. This is the only optical absorption possibility 
for an intrinsic semiconductor NC under equilibrium. By 
relaxing these two constraints we can introduce other ab- 
sorption channels, mainly through carrier injection or op- 
tical pumping. The associated absorption in either case 
is sometimes referred to as "free" carrier absorption de- 
spite the carrier confinement in NCs. In our work we 



discriminate between the two. The electrical injection or 
doping gives rise to intraband absorption, also termed 
as intersubband absorptions^ which has practical impor- 
tance for mid- and near-infrared photodetectors,^^ The 
optical pumping which is usually well above the effective 
band gap leads to excited-state absorption (also termed 
as photoinduced absorption) which is an undesired effect 
that can inhibit the development of optical gain^i Recent 
experiments on excited-state absorption concluded that 
more attention should be devoted to the role of the exci- 
tation conditions in the quest for the silicon laser 
Therefore, the aim of this work is to provide a compre- 
hensive theoretical account of all of these direct photon 
absorption mechanisms in Si and Ge NCs under various 
size, shape and excitation conditions. This provides a 
complementary track to the existing experimental efforts 
where the size and shape control are currently major ob- 
stacles. 

The absorption coefficient of the semiconductor NCs 
depends on the product of the optical transition os- 
cillator strength and their joint density of states as 
well as to their volume filling factor within the ma- 
trix. Therefore, the essential decision on a theoretical 
study is the sophistication level of the electronic struc- 
ture. The usual trade off between the computational 
cost and accuracy is operational. The constraints on 
the former are quite stringent as a NC including the 
active region of the matrix surrounding itself can con- 
tain on the order of ten thousand atoms. As for the 
latter, not only the accuracy but also the validity of a 
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chosen approach can become questionable. Computa- 
tionaUy low-cost approaches like the envelope function in 
conjunction with 8-band k-p are not as accurate for this 
task and furthermore, they miss some critical symme- 
tries of the underlying lattice.'^*' On the other extreme, 
there lies the density functional theory-based ab initio 
codes'^''^ which have been applied to smaller NCs con- 
taining less than 1000 atoms which still require very de- 
manding computational resources i^'^^'^^i^^'^''i^^'^^i'^^ The 
ab initio analysis of larger NCs of sizes between 3- 
10 nm is practically not possible with the current com- 
puter power. While this technological hurdle will be 
gradually overcome in the years to come, there exists 
other atomistic approaches that can be employed for 
NC research which can be run on modest platforms and 
are much simpler to develop, such as the tight binding 
technique which has been successfully employed by sev- 
eral groups.— 1^'^'^ On the pseudopotential-based ap- 
proaches, two new recipes were proposed by Wang and 
Zunger over the last decadei^^i^i^^ The folded spectrum 
method^ relies on standard plane wave basis and di- 
rect diagonalization; its speed is granted from being fo- 
cused on relatively few targetted states. For the study 
of excitons this approach becomes very suitable whereas 
for the optical absorption spectra where a large number 
of states contribute it loses its advantage. Their other 
recipe is the so-called linear combination of bulk bands 
(LCBB)i^i22i^ As a matter of fact, the idea of using 
bulk Bloch states in confined systems goes back to ear- 
lier times, one of its first implementations being the stud- 
ies of Ninno et al^^ Up to now, it has been used for 
self- assembled quantum dots^i^ superlattices^i^ and 
high-electron mobility transistors,— and very recently on 
the nc-Si aggregation stages.— In this work, we apply 
LCBB to the electronic structure and absorption spectra 
of Si and Ge NCs. An important feature of this work, in 
contrast to commonly studied hydrogen-passivated NCs 
is that we consider NCs embedded in a wide band-gap 
matrix which is usually silica4^ In principle, other matri- 
ces such as alumina or silicon nitride can be investigated 
along the same lines. 

The organization of the paper is as follows: in Sec- 
tion II we describe the theoretical framework which in- 
cludes some brief information on the LCBB technique 
and the absorption expressions. A self-critique of the the- 
oretical model is done in Section III. Section IV presents 
the results and discussions on the band edge electronic 
structure, interband, intraband, and excited-state ab- 
sorptions followed by our conclusions in Section V. Ap- 
pendix section contains technical details on the employed 
pseudopotential form factors and our LCBB implemen- 
tation. 



II. THEORY 



method which is particularly convenient for embedded 
NCs containing several thousand atoms.— The fact 
that it is a pseudopotential-based method makes it more 
preferable over the empirical tight binding technique for 
the study of optical properties as aimed in this work. In 
this technique the NC wavefunction with a state label j 
is expanded in terms of the bulk Bloch bands of the con- 
stituent core and/or embedding medium (matrix) mate- 
rials 

n,k.{7 

where N is the number of primitive cells within the com- 
putational supercell, - is the expansion coefficient 

set to be determined and a is the constituent bulk mate- 
rial label pointing to the NC core or embedding medium. 
u'^ -(r) is the cell-periodic part of the Bloch states which 
can be expanded in terms of the reciprocal lattice vectors 
0} as 

G 

where f^o is the volume of the primitive cell. The atom- 
istic Hamiltonian for the system is given by 

where W^(Rj) is the weight function that takes values 
or 1 depending on the type of atom at the position 
Rj — d'^r^ and is the screened spherical pseudopoten- 
tial of atom a of the material a. We use semiempirical 
pseudopotentials for Si and Ge developed particularly for 
strained Si/Ge superlattices which reproduces a large va- 
riety of measured physical data such as bulk band struc- 
tures, deformation potentials, electron-phonon matrix el- 
ements, and heterostructure valence band offsets^i With 
such a choice, this approach benefits from the empirical 
pseudopotential method (EPM), which in addition to its 
simplicity has another advantage over the more accurate 
density functional ab initio techniques that run into well- 
known band-gap problem'^- which is a disadvantage for 
the correct prediction of the excitation energies. 

The formulation can be cast into the following gener- 
alized eigenvalue equationi^Sii^ 

H ,p , t C r = E^ S ,p , ,- C r: (4) 

n,k,(T n^k^a 

where 

Hn'k'a',nka = V |f + Ktal l^fccr) , 



For the electronic structure of large-scale atomistic 
systems Wang and Zunger have developed the LCBB 
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G,G 



G + k-G' -k' 



S .r, , r = (n'k'a'lnka 

n'k' tj' .7ika \ I 



Here, the atoms are on regular sites of the underlying 



Bravais lattice: R„ 



niSi + 712(12 + "-3^3 where 



{Si} are its direct lattice vectors of the Bravais lat- 
tice. Both the NC and the host matrix are assumed to 
possess the same lattice constant and the whole struc- 
ture is within a supercell which imposes the periodic- 
ity condition W (^Rm ,n2 ,n3 

recalling its Fourier representation W (^Rni,n2.na 



implies e 



iq- Ni Si 



1, so that q 



gmi,m2,m3 = + &2 ^ + &3 , whcrc {h} are the 

reciprocal lattice vectors of the bulk material. Thus the 
reciprocal space of the supercell arrangement is not a 
continuum but is of the grid form composed of points 

{(fmi.ma.ma}) whcrC nii ^ 0, I, . . . , Ni - 1. 

An important issue is the choice of the host matrix 
material. If the NC is surrounded by vacuum, this corre- 
sponds to the free-standing case. However, the dangling 
bonds of the surface NC atoms lead to quite a large num- 
ber of interface states which adversely contaminate espe- 
cially the effective band-gap region of the NC. In practice 
NCs are embedded into a wide band-gap host matrix 
which is usually silica4S However, the pseudopotential 
for oxygen is nontrivial in the case of EPM^ and further- 
more, lattice constant of Si02 is not matched to either of 
the core materials introducing strain effects. Therefore, 
we embed the Si and Ge NCs into an artificial wide band- 
gap medium which for the former reproduces the proper 
band alignment of the Si/Si02 system. To circumvent the 
strain effects which are indeed present in the actual sam- 
ples, we set the lattice constant and crystal structure of 
the matrix equal to that of the core material. The pseu- 
dopotential form factors of the wide band-gap matrices 
for Si and Ge can easily be produced starting from those 
of the core materials. More details are provided in the 
Appendix section. The resultant bulk band structures 
for Si and Ge and their host wide band-gap matrices are 
shown in Fig. [TJ With the use of such a lattice-matched 
matrix providing the perfect termination of the surface 
bonds of the NC core atoms lead to the removal of all gap 
states as can be observed in Fig. [5] In these plots, the 
evolution of the effective band-gaps towards their bulk 
values (marked by dashed lines) is clearly seen as the 
diameter increases. 

Once the electronic wavefunctions of the NCs are avail- 
able, their linear optical properties can be readily com- 
puted. The three different types of direct (zero-phonon) 
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FIG. 1: (Color online) EPM band structures for bulk (a) 
Si, (b) Ge together with their wide band-gap matrices (thick 
lines) which for the former reproduces the band line-up of the 
Si/Si02 interface. 



photon absorption processes considered in this work are 
illustrated in Fig. [3] These are interband, intraband and 
excited-state absorptions. In the latter, the blue (dark- 
colored) arrow represents optical pumping and follow- 
ing carrier relaxation, the downward green arrow corre- 
sponds to luminescence which can be to a final interface 
state (dashed line)i^ For all these processes, the relevant 
quantity is the imaginary part of the dielectric function. 
Within the independent-particle approximation and the 
artificial supercell frameworkii it becomes 



Im{eaa('^)} = 



J c 



r/(27r) 



[E, - K„ - H" + (r/2) 



2' 



(5) 



where, a = x,y,z denotes the cartesian components of 
the dielectric tensor and 
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Ec — Ey 



(6) 



is the oscillator strength of the transition. In these ex- 
pressions mo is the free electron mass, e is the magnitude 
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FIG. 2: (Color online) The variation of NC states with respect 
to diameter for Si and Ge NCs. The bulk band edges are 
marked with a dashed line for comparison. 
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FIG. 3: (Color online) Illustration for the three different ab- 
sorption processes in NCs considered in this work; interband, 
intraband and excited-state absorption. The yellow (light- 
colored) arrows indicate the direct photon absorption transi- 
tions, the blue (dark-colored) arrow represents optical pump- 
ing and the downward green arrow corresponds to lumines- 
cence which can be to a interface state (dashed line). 



of the electronic charge, and T is the full-width at half 
maximum value of the Lorentzian broadening. The label 
V (c) correspond to occupied (empty) valence (conduc- 
tion) states referring only to their orbital parts in the 
absence of spin-orbit coupling; the spin summation term 
is already accounted in the prefactor of Eq. [51 Finally, 



Vsc is the volume of the supercell which is a fixed value 
chosen conveniently large to accommodate the NCs of 
varying diameters, however, if one uses instead, that of 
the NC, l^NC, this corresponds calculating Im{eaa}//^ 
where /„ — Vnc/^sc is the volume filling ratio of the 
NC. For the sake of generality, this is the form we shall 
be presenting our results. The electromagnetic intensity 
absorption coefficient a{uj) is related to the imaginary 
part of the dielectric function through^ 



lm{eaa{i^)} = — aaa(w) , 



(7) 



where is the index of refraction and c is the speed of 
light. 

In the case of intraband absorption, its rate depends 
on the amount of excited carriers. Therefore, we consider 
the absorption rate for one excited electron or hole that 
lies at an initial state i with energy Ei. As there are a 
number of closely spaced such states, we perform a Boltz- 
mann averaging over these states as e~^^^ / e~^^i . 
We further assume that the final states have no occu- 
pancy restriction, which can easily be relaxed if needed. 
The expression for absorption rate per an excited carrier 
in each NC becomes 



an 



-l3Ei 



r/(2^) 



[Ef -E,- H + (r/2)2 



(8) 



where again a is the light polarization direction. 

Finally, we include the surface polarization effects, also 
called local field effects (LFE) using a simple semiclassi- 
cal model which agrees remarkably well with more rigor- 
ous treatments.^*' We give a brief description of its im- 
plementation. First, using the expression 



esc — fv^NC + (1 ^ /i))einatrix , 



(9) 



we extract (i.e., de-embed) the size-dependent NC di- 
electric function, enci where esc corresponds to Eq. \S\ 
suppressing the cartesian indices. Cmatrix is the dielec- 
tric function of the host matrix; for simplicity we set 
it to the permittivity value of Si02, i.e., ematrix — 4. 
Since the wide band-gap matrix introduces no absorp- 
tion up to an energy of about 9 eV, we can approxi- 
mate ImjeNc} = Iiii{esc}//i)- One can similarly ob- 
tain the Re{eNc} within the random-phase approxima- 
tion's, hence get the full complex dielectric function eNC • 
According to the classical Clausius-Mossotti approach, 
which is shown to work also for NCs,'^^ the dielectric 
fimction of the NC is modified as 



ENCLFE — Ematrix 



4eNC — Ematrix 
CNC + 2ematrix, 



(10) 



to account for LFE. The corresponding supercell dielec- 
tric function, esc, LFE follows using Eq. [9l Similarly, the 
intensity absorption coefficients are also modified due to 
surface polarization effects, cf. Eq. [T] Its consequences 
will be reported in Section IV. 
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III. A SELF-CRITIQUE OF THE 
THEORETICAL MODEL 



The most crucial simplification of our model is the fact 
that strain-related effects are avoided, a route which is 
shared by other theoretical works JiiSiiS^iS^i^ For large 
NCs this may not be critical, however, for very small 
sizes this simplification is questionable. An important 
support for our act is that Weissker and coworkers have 
concluded that while there is some shift and possibly a re- 
distribution of oscillator strengths after ionic relaxation, 
the overall appearance of the absorption spectra does not 
change stronglyii^ We should mention that Wang and 
Zunger have offered a recipe for including strain within 
the LCBB framework, however, this is considerably more 
involved. ■^^ Another widespread simplification on Si and 
Ge NCs is the omission of the spin-orbit coupling and 
the nonlocal (angular momentum-dependent) pseudopo- 
tential terms in the electronic structure Hamiltonian. Es- 
pecially the former is not significant for Si which is a light 
atom but it can have a quantitative impact on the valence 
states of Ge NCs; such a treatment is available in Ref . Issl. 

On the dielectric response, there are much more sophis- 
ticated and involved treatments^A whereas ours is equiv- 
alent to the independent particle random phase approxi- 
mation^*^ of the macroscopic dielectric function with the 
surface polarization effects included within the classical 
Clausius-Mossotti modeli^ The contribution of the ex- 
cluded excitonic and other many-body effects beyond the 
mean-field level can be assessed a posteriori by compar- 
ing with the available experimental data. However, it 
is certain that the precedence should be given to classi- 
cal electrostatics for properly describing the background 
dielectric mismatch between the core and the wide band- 
gap matrix.^*' In our treatment this is implemented at an 
atomistic level. 

Another effect not accounted in this work is the role 
of the interface region. Our wide band-gap matrix can 
reproduce the proper band alignment and dielectric con- 
finement of an Si02 matrix, however, the interface chem- 
istry such as silicon-oxygen double bonds^ are not repre- 
sented. These were shown to be much more effective on 
the emission spectrai^ii Nevertheless, our results can be 
taken as the benchmark for the performance of the atom- 
istic quantum and dielectric confinement with a clean and 
inert interface. Finally, we do not consider the phonon- 
assisted^ or nonlinear absorption. The list of these ma- 
jor simplifications also suggest possible improvements of 
this work. 



IV. RESULTS AND DISCUSSIONS 



standing. However, we first begin with the dependence 
of the optical gap on the NC size, mainly as a check of 
our general framework. There exist two different atomic 
arrangements of a spherical NC depending on whether 
the center of the NC is an atomic position or a tetrahe- 
dral interstitial location; under no ionic relaxation, Del- 
ley and Steigmeier have treated both of these classes as 
having the Td point symmetry^SS However, the tetrahe- 
dral interstitial-centered arrangement should rather have 
the lower point symmetry of C^y and it is the arrange- 
ment that we construct our NCs. This leads to even 
number of NC core atoms, whereas it becomes an odd 
number with the Td point symmetry. We identify the ir- 
reducible representation of a chosen NC state by checking 
its projection to the subspace of each representation!^ 
For the C^v point group these are denoted by Ai, A2, 
and E. We utilize this group-theoretic analysis in the 
next subsections. 
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FIG. 4: (Color online) Comparison of optical gap as a func- 
tion of NC diameter of this work with previous experimental 
and theoretical data: Furukawa,^" Kanemitsu,Ji Takeoda,— 
Ogiit,-^ Vasiliev,-i& Garoufalis,-^ Niquet,-^ Tsolakidis.— 



In this section we present our theoretical investigation 
of the linear optical properties of Si and Ge NCs. Three 
different direct photon absorption processes are consid- 
ered as illustrated in Fig. [3] each of which can serve for 
technological applications as well as to our basic under- 



A. Effective optical gap 

The hallmark of quantum size effect in NCs has been 
the effective optical gap with quite a number of theoret- 
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studies per- 



formed within the last decade. Figure [J] contains a com- 
pilation of some representative results. For Si NCs, it can 
be observed that there is a good agreement among the 
existing data, including ours. On the other hand, for the 
case of Ge NCs there is a large spread between the exper- 
imental data whereas our theoretical results are in very 
good agreement with both ab initic^'^ and tight binding 
resultsi^ In our approach the optical gap directly corre- 
sponds to the LUMO-HOMO energy difference, as calcu- 
lated by the single-particle Hamiltonian in Eq. ([3]). This 
simplicity relies on the finding of Delerue and coworkers 
that the self-energy and Coulomb corrections almost ex- 
actly cancel each other for Si NCs larger than a diameter 
of 1.2 nm.^ 
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FIG. 5: (Color online) The variation of HOMO and LUMO 
energies with respect to NC diameter for Si and Ge NCs that 
belong to Czv point group. 
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FIG. 6; (Color online) (a) The energy difference between dou- 
bly degenerate and nondegenerate states, one of which be- 
comes the HOMO with respect to diameter in Si and Ge NCs 
that belong to Czv point group; solid lines are for guiding the 
eyes; ao is the lattice constant for Si or Ge NC. (b) The low- 
est three conduction states, one of which becomes the LUMO 
with respect to diameter in Si NCs. 



B. HOMO and LUMO oscillations with respect to 

size 

When we plot the variation of individual LUMO and 
HOMO levels as in Fig.[5]we observe with the exception of 
nc-Ge LUMO curve some non-smooth behavior that gets 
pronounced towards smaller sizes. The triple degeneracy 
in the absence of spin-orbit coupling of the valence band 
maximum in bulk Si and Ge is lifted into two degener- 
ate and one nondegenerate states. The energy difference 
between these two set of states is observed to display 
an oscillatory behavior as the NC size gets smaller as 
shown in Fig. [S] (a). Using the C31, point group symme- 
try operations we identify the doubly degenerate states 
to belong to E representation and nondegenerate state to 



A\ or Ai- Furthermore, we observe a similar oscillation 
in the LUMO region of Si NCs as shown in Fig. [6l^b). 
The low-lying conduction states of Si NCs form six-pack 
groups which is inhereted from the six equivalent 0.85X 
conduction band minima of hulk Si. The confinement 
marginally lifts the degeneracy by sampling contributions 
from other parts of the Brillouin zone. This trend is ob- 
served in Fig. mb) as the NC size gets smaller. On the 
other hand, for Ge NCs all LUMO states belong to the 
same A\ representation and therefore shows no oscilla- 
tions (cf. Fig. [5]). Ultimately, the source of these oscilla- 
tions is the variation of the asphericity of the NCs of C3.U 
point symmetry with respect to size, which can energeti- 
cally favor one of the closely spaced states. In the case of 
the LUMO state of Ge NCs, there is a substantial energy 
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gap between LUMO and the next higher-lying state. 

For further insight, we display in Fig. [7] the isosur- 
face plots of the envelope of the six highest states up to 
HOMO for a Si NC of diameter 2.16 nm. Point group 
representation of each state is also indicated. For this 
particular diameter, HOMO has E representation which 
is twofold degenerate. The nondegenerate Ai state also 
becomes the HOMO for different diameters. This is illus- 
trated in Table I which shows the evolution of the HOMO 
and LUMO symmetries as a function of diameter for Si 
and Ge NCs. There, it can be observed that for the latter 
the HOMO can also acquire the A2 for larger diameters. 



that its size dependence is much stronger than the energy 
dependence. This reduction in the absorption due to LFE 
can become a major concern for solar cell applications. 
It needs to be mentioned that this effect is highly sensi- 
tive to the permittivity mismatch between the core and 
matrix media. To illustrate this point, in Fig. [5] the case 
for AI2O3 matrix (having a permittivity of 9.1) is also 
displayed for 1.41 nm Si NC, where it can be seen that 
compared to Si02 (with a permittivity of about 4) the 
reduction in absorption due to LFE is much less. Based 
on this finding, we employ these size-dependent absorp- 
tion reduction factors in the results to follow including 
the intraband and excited-state cases. 
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FIG. 7: (Color online) The isosurfaces of the envelopes of the 
wavefunctions of the highest-lying six valence states up to 
HOMO for a Si NC of diameter 2.16 nm. The isosurfaces are 
drawn for the 95% of the peak value of the envelope wavefunc- 
tions. The Csv point group representations are indicated for 
each wavefunction. Note that some of the plots are rotated 
with respect to others for best viewing angle. 



C. Interband absorption 

The interband absorptions of Si and Ge NCs for a va- 
riety of diameters are shown in Fig. [S] For a fair com- 
parison, all different size NCs should possess the same 
volume filling factor. Therefore, we display the results at 
unity volume filling or equivalently per ^ The left and 
right panels display the cases without and with surface 
polarization effects (or LFE), respectively. There exists 
remarkable differences between the two for both Si and 
Ge NCs. For instance, even though Ge NCs do not show 
significant size dependence without LFE, this is not the 
case when LFE is included. From the ratio of both pan- 
els, the so-called local field absorption reduction factor 
can be extracted as shown in Fig. [51 It can be observed 




Photon Energy (eV) 
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FIG. 8: (Color online) The imaginary part of dielectric func- 
tion for unity volume filling factor, /„ for Si and Ge NCs at 
different diameters with (right panel) and without (left panel) 
local field effects. A Lorentzian broadening energy full width 
of 200 meV is used. 



In Fig. [To] we compare our results with the experimen- 
tal data of Wilcoxon et al. for Si NCsH and Ge NCs 
There is a good overall agreement in both cases espe- 
cially with LFE, however, for the case of Si NCs this is 
much more satisfactory. The major discrepancies can be 
attributed to excitonic effects that are not included in 
our work. In the case of Si NCs (Fig. fTUTa)). we also 
display the tight binding result of Trani et al. which also 
includes LFE.^" 

An issue of practical concern is the effect of deviation 
from the spherical shape of the NCs depending on the 
growth conditions. At this point we would like to in- 
vestigate the effect of shape anisotropy on the interband 
absorption. Starting from a spherical 2.54 nm diameter 
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TABLE I: Czv irreducible representations of the HOMO and LUMO in Si and Ge NCs of various diameters {D). 
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D (nm) 


HOMO 
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Ai 


Ai 


38 


1.13 


Ai 


Ai 


1.18 


Ai 


Ai 


56 


1.29 


Ai 


Ai 


1.34 


Ai 


Ai 


74 


1.41 
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FIG. 9; (Color online) The local field absorption factor ex- 
tracted from the previous figure. 



Ge NC, we form prolate and oblate ellipsoidal NCs with 
ellipticities e = —0.6 and +0.6, respectively. All three 
NCs contain the same number of 384 core atoms; the 
atomic arrangement of the ellipsoidal NCs are shown in 
the inset of Fig.[TT] In the same figure we compare the zz 
components of the imaginary part of the dielectric ten- 
sor for three different ellipticities. It is observed that the 
effect on the interband absorption is not significant; the 
difference is even less for the Si NCs (not shown). 



D. Intraband absorption 

Unlike the interband case, for the intraband absorption 
we need to introduce electrons to the conduction states 
or holes to the valence states by an injection mechanism. 
We assume that after injection these carriers relax to 
their respective band edges and attain a thermal distri- 
bution. Therefore, we perform a Boltzmann averaging at 
room temperature (300 K) over the initial states around 
LUMO (HOMO) for electrons (holes). The absorption 
coefficients to be presented are for unity volume filling 
factors and for one carrier per NC; they can easily be 
scaled to a different average number of injected carriers 
and volume filling factors. In Fig.[T2]the Si NCs of differ- 
ent diameters are compared. The intraband absorption 
is observed to be enhanced as the NC size grows up to 
about 3 nm followed by a drastic fall for larger sizes. For 
both holes and electrons very large number of absorp- 
tion peaks are observed from 0.5 eV to 2 eV. Recently, 
de Sousa et al. have also considered the intraband ab- 
sorption in Si NCs using the effective mass approximation 
and taking into account the multi- valley anisotropic band 
structure of Si.^^ However, their absorption spectra lacks 
much of the features seen in Fig.[T3 Turning to Ge NCs, 
shown in Fig. [13] the intravalence band absorption profile 
is very similar to that of Si NCs, however in this case the 
intraconduction band absorption is much weaker. 

Mimura et al. have measured the optical absorption 
in heavily phosphorus doped Si NCs of a diameter of 
4.7 nm.®^ This provides us an opportunity to compare our 
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FIG. 10: (Color online) Comparison of our absorbance re- 
sults with the available data: for Si, experimental work of 
Wilcoxon^"^ and the theoretical tight binding results of Tran;"'" 
and for Ge, the experimental work of Wilcoxon^^. For our 
spectra a Lorentzian broadening energy full width of 200 meV 
is used. 

results on the intraconduction band absorption in Si NCs. 
There is a good order-of-magnitude agreement. However, 
in contrast to our spectra in Fig. 1 121 which contains well- 
resolved peaks, they have registered a smooth spectrum 
which has been attributed by the authors to the smearing 
out due to size and shape distribution within their NC 
ensemble;^ 



E. Excited-state absorption 

Finally, we consider another intraband absorption pro- 
cess where the system is under a continuous interband 



tz nc-Ge 
D=2.54 nm 




Photon Energy (eV) 



FIG. 11: (Color online) The effect of ellipticity, e on the 
lm{ezz} / fv for a Ge NC with a diameter of 2.54 nm. The 
insets show NC core atoms of the prolate (e = —0.6) and 
oblate (e = 0.6) shapes; the z-direction is also indicated. A 
Lorentzian broadening energy full width of 200 meV is used. 

optical pumping that creates electrons and holes with ex- 
cess energy. We consider three different excitation wave- 
lengths: 532 nm, 355 nm and 266 nm which respectively 
correspond to the second-, third- and fourth- harmonic of 
the Nd-YAG laser at 1064 nm. The initial states of the 
carriers after optical pumping are chosen to be at the pair 
of states with the maximum oscillator strength*** among 
interband transitions under the chosen excitation. The 
determined energies of these states are tabulated in Ta- 
ble II where it can be observed that in general the excess 
energy is unevenly partitioned, mainly in favor of the 
conduction states. Once again a Boltzmann averaging is 
used to get the contribution of states within the thermal 
energy neighborhood. 

Considering 3 nm and 4 nm diameters, the results are 
shown in Figs . fT4l and [TSl for Si and Ge NCs, respectively. 
Note that the 532 nm excitation results are qualitatively 
similar to those in intraband absorption, cf. Figs. [T^] 
and [ini This is expected on the grounds of small ex- 
cess energy for this case. Some general trends can be 
extracted from these results. First of all, the conduction 
band absorption is in general smooth over a wide en- 
ergy range. On the other hand the valence band absorp- 
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FIG. 12: (Color online) Intravalence and intraconduction 
state absorption coefficients in Si NCs of different diameters 
per excited carrier and at unity filling factor. A Lorentzian 
broadening energy full width of 30 meV is used. Mind the 
change in the vertical scale for 4 nm diameter case. 



tion contains pronounced absorption at several narrow 
energy windows mainly below 1 eV and they get much 
weaker than the conduction band absorption in the re- 
maining energies. As the excitation energy increases the 
absorption coefRcient per excited carrier in general de- 
creases. In connection to silicon photonics, we should 
point out that the excited-state absorption is substantial 
including the important 1.55 /um fiber optics communi- 
cation wavelength. These results provide a more com- 
prehensive picture than the reported experimental mea- 
surementa^i^^i^S which are usually obtained at a single 
energy of the probe beam. Finally, it needs to be men- 
tioned that for both intraband and excited-state absorp- 
tions displayed in Figs. [12] to [151 the high energy parts 
will be masked by the interband transition whenever it 
becomes energetically possible. 



V. CONCLUSIONS 

The subject of Si and Ge NCs has become an estab- 
lished research field. A fundamental process such the 



FIG. 13: (Color online) Same as Fig. [12] but for Ge NCs. 



direct photon absorption deserves further investigation 
from a number of perspectives. In this theoretical study, 
we consider the interband, intraband and excited-state 
absorption in embedded Si and Ge NCs of various sizes. 
For this purpose, we developed an atomistic pseudopo- 
tential electronic structure tool, the results of which agree 
very well with the published data. It is further observed 
that the HOMO for both Si and Ge NCs and the LUMO 
for Si NCs display oscillations with respect to size among 
different representations of the C31, point group to which 
these spherical NCs belong. Our detailed investigation 
of the interband absorption reveals the importance of 
surface polarization effects that significantly reduce the 
absorption when included. This reduction is found to 
increase with decreasing NC size or with increasing per- 
mittivity mismatch between the NC core and the host 
matrix. These findings should be taken into account 
for applications where the absorption is desired to be 
either enhanced or reduced. For both NC types the de- 
viation from sphericity shows no pronounced effect on 
the interband absorption. Next, the intraband process is 
considered which has potential applications on mid- and 
near-infrared photodetection. The intravalence band ab- 
sorption is stronger compared to intraconduction band 
especially below 1 eV. Finally, we demonstrate that op- 
tical pumping introduces a new degrees of freedom to 
intraband absorption. This is studied under the title of 
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TABLE II: The excited-state energies of the carriers within the valence and conduction states under three different interband 
pump energies for nc-Si and nc-Ge. The energies are given in eV and measured from the HOMO and LUMO, respectively. 
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FIG. 14: (Color online) Excited-state absorption within va- 
lence and conduction states of Si NCs per excited carrier and 
at unity filling factor under three different optical pumping 
wavelengths of 532 nm, 355 nm and 266 nm. Dotted lines 
in black color refer to total absorption coefficients. Two dif- 
ferent diameters are considered, 3 and 4 nm. A Lorentzian 
broadening energy full width of 30 meV is used. 



FIG. 15: (Color online) Same as Fig.[T4]but for Ge NCs. 
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TABLE III: Parameters of the pseudopotential form factors of Si, Ge and their wide band-gap matrices, ao is the lattice 
constant. See text for the units. 





ao (A) 


ai 


a.2 




04 


as 


oe 


Si 


5.43 


1.5708 


2.2278 


0.606 


-1.972 


5.0 


0.3 


Matrix-Si 


5.43 


1.5708 


2.5 


0.135 


-13.2 


6.0 


0.3 


Ge 


5.65 


0.7158 


2.3592 


0.74 


-0.38 


5.0 


0.3 


Matrix-Ge 


5.65 


0.4101 


2.7 


0.07 


-2.2 


5.0 


0.3 



APPENDIX: SOME TECHNICAL DETAILS ON 
THE LCBB IMPLEMENTATION 

In this section we first provide the details on the pseu- 
dopotential form factors of the bulk Si, Ge and their as- 
sociated wide band-gap matrices to substitute for Si02. 
We use the local empirical pseudopotentials for Si and 
Ge developed by Friedel, Hybertsen and Schliiter.— They 
use the following functional form for the pseudopotential 
form factor at a general wave number q: 

y.Al) - W tanh ( + i" . 

^^^^^ e°3(9'-"4) + 1 [2 V "6 / 2_ 

Using the parameters supplied in Table III, the pseudopo- 
tential form factors come out in Rydbergs and the wave 
number in the above equation should be taken in atomic 
units (1/Bohr radius). Another important technical re- 
mark is about the EPM cut off energies. We observed 



that even though the EPM band energies (i.e., eigen- 
values) converge reasonably well with cut off energies as 
low as 5-10 Ry, the corresponding Bloch functions (i.e., 
eigenvectors) require substantially higher values to con- 
verge,— The results in this study are obtained using 14 
and 16 Ry for Si, Ge, respectively. 

Finally, some comments on the LCBB basis set con- 
struction is in order. We only employ the bulk bands of 
the core material. The bulk band indices are chosen from 
the four valence bands an the lowest three to four con- 
duction bands; usually these are not used in conjunction 
but separately for the NC valence and conduction states, 
respectively. The basis set is formed from a sampling 
over a three-dimensional rectangular grid in the recip- 
rocal space centered around the F-point. Its extend is 
determined by the full coverage of the significant band 
extrema, such as conduction band minima of Si at the 
six equivalent 0.85 X points. The final LCBB basis set 
typically contains some ten thousand members. 
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